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1. Introduction 

Ultra-relativistic heavy ion collisions offer the unique opportunity to probe highly excited 
dense nuclear matter under controlled laboratory conditions. One of the main driving 
forces for these studies is the expectation that an entirely new form of matter may be 
created in such reactions. This form of matter, called the Quark Gluon Plasma (QGP), is 
the QCD analogue of the plasma phase of ordinary atomic matter However, unlike such 
ordinary plasmas, the deconfined quanta of a QGP are not directly observable because of 
the fundamental confining property of the physical QCD vacuum. What is observable are 
hadronic and leptonic residues of the transient QGP state. The QGP state formed in nu- 
clear collisions is a transient rearrangement of the correlations among quarks and gluons 
contained in the incident baryons into a larger but globally still color neutral system with 
however remarkable theoretical properties, such as restored chiral symmetry and decon- 
finement. The task with heavy ion reactions is to provide experimental information on 
this fundamental prediction of the Standard Model (for recent reviews on QGP signatures, 
please see [1]). 



2. The QCD Phase Diagram 

QCD is a non-abelian gauge theory, it's basic constituents are quarks and anti-quarks in- 
teracting through the exchange of color-charged gluons. At very high temperatures and 
densities, in the domain of weak coupling between quarks and gluons, long range inter- 
actions are dynamically screened [2,3]. Quarks and gluons are then no longer confined to 
bound hadronic states ("deconfinement"). Furthermore, chiral symmetry is restored - for 
baryon-free matter - apparently at the same temperature Tc- 

Lattice gauge simulations of QCD [4,5], provide the only rigorous method to compute 
the equation of state of strongly interacting elementary particle matter In principle both, 
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3-flavour phase diagram 



|t;f" ~ 175 MeV| ^^^^ p^,^ |Tj ~ 270 MeV| 




Figure 1. Left frame: QCD phase diagram of 3-flavor QCD witii degenerate 
(u,d)-quark masses and a strange quark mass rris. Right frame: schematic QCD phase 
diagram as a function of /ib- 



the non-perturbative hadronic matter and the non-perturbative QGP phases of QCD can be 
investigated. The main disadvantage of lattice simulations is the current practical restric- 
tion to finite, periodic, baryon free systems in global equilibrium, a scenario far away from 
the highly inhomogeneous off-equilibrium situation found in complex heavy-ion reactions. 

Lattice calculations yield a critical temperature of Tc « 270 MeV in the quenched 
approximation [6,7] - where neither dynamical quarks, nor a chiral phase transition exist. 
Simulations including dynamical quarks at /i^ =0 indicate a critical temperature in the 
order of Tc ~ 165 MeV (see figure 1, taken from [7]). 

Many striking QGP signatures depend heavily on the assumption of a first order phase 
transition and the existence of a mixed phase of QCD matter However, both, the order 
of the phase transition as well as the critical temperature, depend on the parameters of the 
calculation, namely the number of quark flavors and their masses. For the most realistic 
case of QCD with two flavors of light quarks with masses between 5 and 10 MeV and 
one flavor with a mass around 200 MeV, the situation remains unclear: the order of the 
phase transition seems to depend on the numerical values for the masses of the light and 
heavy quarks [8,7]. If the latter is too heavy, the transition might be smeared out to a mere 
rapid increase of the energy density over a small temperature interval. In this case the use 
of simple deconfinement scenarios may lead to wrong expectations for observables. The 
elementary excitations in such a phase transition region ought not be described by quarks 
and gluons but could physically resemble more hadronic excitations with strongly modified 
"in-medium" properties [9]. 

This raises a practical question, whether conclusions based on /i^ = estimates, might 
misguide physical argumentation for observables in nuclear collisions. This warning is 
particularly appropriate for those QGP-signals, where a 50% quantitative change of an ob- 
servable is used to differentiate QGP production scenarios from ordinary hadronic transport 
ones. 

The inclusion of the second most important thermodynamic variable, the chemical po- 
tential fiB into a full fletched lattice-QCD calculation has recently been attempted [10], but 
still carries large systematical uncertainties. These calculations predict the existence of a 
tricritical point - only for fis > I^B.crit would a phase-boundary with a possibly 1st order 
phase-transition exist. 

An additional complication is that for systems of finite volume (V< 125 fm'^) the de- 
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Figure 2. Transport theory approaches for RHIC and their range of apphcabihty. 



confinement cannot be complete. Fluctuations lead to a finite probability of the hadronic 
phase above Tc- The sharp discontinuity (e.g. e/T'^) is thus smeared out [11]. 

Recently, a lot of progress has been made towards understanding the QCD phase di- 
agram at large chemical potential [12]. At very high densities, asymptotic freedom 
and BCS theory have been utilized to predict new, superconducting, phases of QCD mat- 
ter. Depending on the mass of the strange quark (m^ 3> m^^), first a two flavor color- 
superconductor {rus oo) or a color-flavor locked superconducting phase (m^ decreas- 
ing with fiB) are found (see right frame of figure 1, taken from [13]). The color-flavor 
locking refers to the symmetry of the ground state only being given under a simultaneous 
transformation of color and flavor. It should be noted, however, that these new phases of 
high-density QCD will most likely not be found in heavy-ion reactions (the temperature 
in such collisions is too high for color superconductivity to occur), but may rather have 
a large impact in our understanding of matter in neutron stars and other cold super-dense 
objects (for a review on the topic, please see [13]). 

3. Dynamics: Transport Theory at SPS and RHIC 

Transport Theory offers the unique capability of connecting experimentally observable 
quantities in a relativistic heavy ion collision with its time evolution and reaction dynamics, 
thus giving crucial insights into the possible formation of a transient deconfined phase of 
hot and dense matter, the Quark-Gluon-Plasma. 

Figure 2 provides an overview of different transport theoretical ansatzes currently in the 
literature for the description of a relativistic heavy ion collision at RHIC energies. The 
timeline shows a best case scenario for what to expect: the formation of a thermalized 
QGP with subsequent hadronization and freeze-out. Bands with solid lines denote the safe 
range of applicability for the respective transport approach, whereas dashed/dotted bands 
refer to areas in which the approach is still applied but where the assumptions on which 
the approach is based upon may be questionable or not vaUd anymore. 
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The initial state and early pre-equilibrium phase are best described in Classical Yang- 
Mills theory (CYM) [14J or Lattice Gauge Transport (LGT) [15] calculations - only these 
classes of models treat the coherence of the initial state correctly, but do not provide any 
meaningful dynamics for the later reaction stages. 

The Parton Cascade Model (PCM) [16] and related pQCD approaches [33] treat the ini- 
tial state as incoherent parton configuration, but are very well suited for the pre-equilibrium 
phase and subsequent thermalization, leading to a QGP and hydrodynamic expansion. Mi- 
croscopic degrees of freedom in this ansatz are quarks and gluons which are propagated 
according to a Boltzmann Equation with a collision term using leading order pQCD cross 
sections. Augmented with a cluster hadronization ansatz the PCM is applicable up to 
hadronization. The range of this model can be even further extended if it is combined with 
a hadronic cascade which treats the hadronic phase and freeze-out correctly [18]. 

Nuclear Fluid Dynamics (NFD, see e.g. [19,20]) is ideally suited for the QGP and hy- 
drodynamic expansion reaction phase, but breaks down in the later, dilute, stages of the 
reaction when the mean free paths of the hadrons become large and flavor degrees of 
freedom are important. The reach of NFD can also be extended by combining it with a 
microscopic hadronic cascade model - this kind of hybrid approach (dubbed hydro plus 
micro) was pioneered in [21] and has been now also taken up by other groups [22]. It is 
to date the most successful approach for describing the soft physics at RHIC. The biggest 
advantage of NFD is that it directly incorporates an equation of state as input - one of it's 
largest disadvantages is that it requires thermahzed initial conditions and one is not able to 
do an ab-initio calculation. 

Last but not least, string and hadronic transport models [23,24] have been very suc- 
cessful in the AGS and SPS energy domains. They treat the early reaction phase as a 
superposition of hadron-hadron strings and are thus ill suited to describe the microscopic 
reaction dynamics of deconfined degrees of freedom. In the later reaction stages, however, 
they are the best suited approach, since they incorporate the full spectrum of degrees of 
freedom of a hadron gas (including flavor dependent cross sections). 

However, it is important to note that there is not a single transport theoretical ansatz 
currently available, which is able to cover the entire time-evolution of a colUsion at RHIC 
in one self-consistent approach. 

4. SPS: highlights and open questions 

During the last eight years the heavy ion research at the CERN-SPS has succeeded to 
achieve the measurement of a wide spectrum of observables relevant for a QGP search. 
Table 1 provides an overview of these measurements together with an assessment if the 
results would be compatible with the possible creation of a new form of deconfined matter 
[25,26]. In the following sections, several of these results wiU be highlighted in greater 
detail: 



4.1 Multi-strange Baryons 

The relative enhancement of strange and especially multistrange baryons with respect to 
peripheral (or proton induced) interactions has been suggested as a signature for the tran- 
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Table 1. Overview of concepts relevant for a QGP search and their corresponding 
observables investigated at the CERN/SPS. The last column indicates the success of 
the measurement and its compatibility with the possible creation of a new form of 
deconfined matter. 

sient existence of a QGP-phase [27-29]: the main argument being that the (chemical or 
flavor) equilibration times should be much shorter in the plasma phase than in a thermally 
equilibrated hadronic fireball of T ~ 160 MeV. 

The dominant production mechanism in an equilibrated (gluon rich) plasma phase, 
namely the production of ,ss pairs via gluon fusion (gg ss) [27], should allow for equi- 
libration times similar to the interaction time of the colliding nuclei, and to the expected 
plasma lifetime (a few fm/c). 

The yields of strange baryons per event calculated in a microscopic string/hadron trans- 
port model (UrQMD) are shown in the left frame of figure 3 as a function of the number 
of participants for Pb-nPb and pH-Pb collisions at 160 GeV/u [30]. The A + A- (circles), 
S~ + E~- (squares), and + (triangles) values are shown. The stars correspond 
to experimental data of the WA97 collaboration [31]. Open symbols represent the results 
of the standard UrQMD calculations, whereas full symbols exhibit a calculation with an 
enhanced string tension of k = 3GeV/fm, for the most central collisions (A^part > 300). 
Obviously the standard UrQMD calculation, which can be seen as a baseline of known 
hadronic physics, strongly underestimates the (multi-)strange particle yields in central col- 
lisions, in particular in the case of the fl^ . Only the inclusion of non-hadronic medium 
effects, like color-ropes |32], which are simulated by increasing the string-tension for cen- 
tral colUsions, enhances the yield to a level of near-compatibility with the data. Similar 
findings have also been made in the context of HIJING calculations |33|. While these 
findings by no means prove the validity of the color-rope approach, they clearly show the 
necessity of some kind of medium effect beyond regular binary hadronic (re) scattering in 
order to understand the data. Recently, hadronic multi-particle interactions in the early 
dense reaction phase have been suggested to significantly enhance the yield of anti-protons 
and (anti-)hyperons [34,35]. It remains to be seen, however, whether these effects are suf- 
ficient to explain the observed 0~ enhancement or other non-hadronic (i.e. deconfinement 
based) effects need to be taken into account. 
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Figure 3. Left: excitation function of strange baryon multiplicity vs. number of partic- 
ipants in a microscopic transport approach. Right: Hydrodynamical analysis of direct 
photon production at the SPS. 



4.2 Penetrating Probes: Dileptons and Direct Photons 

Direct (tliermal) piiotons in a QGP are created dominantly via qq jg (annihilation) and 
gq 'yq (Compton scattering). The production rate and the momentum distribution of the 
photons depend on the momentum distributions of quarks, anti-quarks and gluons in the 
plasma. If the plasma is in thermodynamic equilibrium, the photons may carry information 
on this thermodynamic state at the moment of their production [36-39]. 

The main hadronic background processes to compete against are pion annihilation tttt 
7/5 and Compton scattering irp 'yn [36,40]. The broad oi resonance may act as an 
intermediate state in scattering and thus provide an important contribution [41,40] via 
it's decay into 'yn. In the vicinity of the critical temperature Tc a hadron gas was shown to 
"shine" as brightly (or even brighter than) a QGP [36]. 

Hydrodynamical calculations can be used to compare purely hadronic scenarios of pho- 
ton radiation with scenarios involving a first/second order phase transition to a QGP. They 
show a reduction in the temperature of the photon spectrum in the event of a first order 
phase-transition [42-44]. The right frame of figure 3 shows a hydrodynamical calculation 
assuming the creation of a QGP and subsequent hadronic rescattering [45]. Recent esti- 
mates of photon production in quark-matter (at two loop level) along with the dominant 
reactions in the hadronic matter leading to photons are used. About half of the radiated 
photons are seen to have a thermal origin. The same treatment and the initial conditions 
provide a very good description to hadronic spectra measured by several groups, lending 
additional support to the conclusion that quark gluon plasma has been formed in these 
collisions. 

Dileptons carry information on the thermodynamic state of the medium at the moment of 
production in the very same manner as the direct photons - since the dileptons interact only 
electromagnetic ally they can leave the hot and dense reaction zone basically undistorted, 
too. It has been found that at the SPS the invariant mass spectrum of low-mass dileptons is 
dominated by pion-pion scattering. Medium modifications to the form-factor are necessary 
for a sohd understanding of the data - the current status of experiment and theory point 
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Figure 4. Left: statistical model fit to hadron ratios measured at the SPS. Right: sta- 
tistical model analysis of the central cell of a heavy-ion collision calculated with a 
microscopic transport model. 

toward a strong broadening of the p resonance without a significant mass-shift toward 
lower masses [46] . 

4.3 Chemical Equilibration 

The success of statistical models in describing the (strange) hadron abundances and ratios 
at the CERN/SPS [47] (see the left frame of figure 4, taken from [48]) and the extracted 
7s values close to 1 have led to the common belief that chemical freeze-out in heavy-ion 
reactions at the SPS occurred very close to - or even at hadrochemical equilibrium and that 
this state most likely has been created by a hadronizing QGP [25,26,49]. 

However, the estimated chemical equilibration times may not be sufficiently rapid to 
cause chemical equilibration before hadronization: calculations based on boost-invariant 
hydrodynamics with rate-equations for quark production [50-52], pQCD rate-equations 
[53] or the Parton Cascade Model [54] all indicate that chemical equilibration (and 
strangeness saturation) cannot be achieved during reaUstic life-times of the deconfined 
phase. It has been suggested (e.g. in [50,51]) that the system would be driven toward and 
come close to chemical equilibrium in the subsequent hadronic phase - a scenario which 
would help to bridge the gap between the calculations indicating insufficient equilibration 
time in the plasma phase and the SPS data apparently close to chemical equilibrium at 
chemical freeze-out. Recent calculations assuming a hadronizing QGP out of chemical 
equilibrium with subsequent hadronic rescatting have shown that rescatting via binary col- 
lisions in the hadronic phase is insufficient to drive the system toward chemical equilibrium 
before the expansion of the system leads to chemical freeze-out [55]. However, hadronic 
multi-particle collisions shortly after hadronization may have a large impact and drive the 
system more effectively toward chemical equilibrium [34,35] - a detailed calculation to 
verify/falsify this speculation has yet to be performed. 

Another question which needs to be raised is whether the temperature and chemical 
potential extracted from the statistical model fits to final state hadron yields or ratios really 



Pramana - J. Phys. 



7 



Stejfen A. Bass 



reflect the thermodynamic state of the system at one particular time during its evolution 
(i.e. the conditions at chemical freeze-out) or are rather the result of a superposition of 
different states, due to individual hadron species decoupling continuously from the system 
(as would be expected from their different mean free paths). The latter view is supported by 
a transport model analysis of the time-evolution of the temperature and chemical potential 
in the central cell of a heavy-ion reaction (see the right frame of figure 4, taken from [56]). 
However, true progress on this issue can only be achieved if a model-independent method 
is found to determine the chemical decouphng time of individual hadron species. 

SPS: final remarks 

The lessons learned from the SPS experiments may be summarized by the following core 
observations ^: 

• SPS experiments have created a new state of high energy-density and temperature 
matter 

• sub-hadronic degrees of freedom offer a viable explanation for many of the observed 
phenomena 

• no single observable or measurement is capable of unambiguous proof for the onset 
of deconfinement 

• only a combination of observables and experiments will be able to deliver proof of a 
new phase of matter 

• the concept of a QGP needs to be rethought: in the SPS domain (and perhaps even at 
RHIC) searching for a fireball of quarks and gluons interacting dominantly according 
to leading order pQCD is unrealistic - most likely deconfinement will occur and 
be observed in the form of a fast evolving phase of partonic degrees of freedom 
interacting in the domain of soft non-perturbative QCD. 

5. RHIC: first results and expectations 

The Relativistic Heavy Ion Collider has provided an impressive array of data in its first year 
of operation and QGP theory faces a stiff challenge to keep up with the pace at which new 
and surprising results emerge. In the following sections a brief overview on some of the 
most exciting and challenging facets of RHIC theory will be given - a more comprehensive 
review on the current status can be found in [57]. 



'^a detailed discussion of charmonium suppression can be found at the end of the RHIC section 
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5.1 Initial Conditions 

An issue of major importance for understanding the physics at RHIC is how to describe the 
initial partonic configuration of the system. In recent years the Color Glass Condensate 
[58] has emerged as one of the leading concepts in that domain: with decreasing x (the 
fraction of the longitudinal momentum of the parent hadron), the gluon density in the 
hadron wavefunction grows faster than the quark density, giving rise to a high-density 
multiparticle state of gluons. Eventually, the increase in the number of gluons has to cede 
and saturation, i.e. a limitation of the maximum gluon density per unit phase-space, occurs. 
Saturation can be characterized by an intrinsic momentum scale Qs (the saturation scale), 
below which non-linear effects in the parton structure function evolution are expected to 
slow down and eventually saturate the increase of the gluon density. The saturation scale 
Qs may also provide a natural cut-off for pQCD based calculations, such as in a parton 
cascade model. As we shall see in the following sections, the Color Glass Condensate 
(CGC) approach provides a viable explanation for a number of observations made at RHIC. 

5.2 Hard Processes 

One of the first questions to ask concerning the initial RHIC results is whether evidence for 
hard processes can be seen and if so, to what extent they dominate particle production and 
the underlying reaction dynamics: the PHOBOS data on the charged particle multiplicity 
at midrapidity and its centrality dependence [59] has been analyzed in a variety of different 
approaches, spanning from saturation model approaches to pQCD noini-jet production via 
hard processes [60-62]. Of particular note is a simple analysis carried out in the eikonal 
approach, decomposing particle production into a hard component, which scales with the 
number of colhsions {Ncoii), and a soft component, which scales with the number of par- 
ticipants {Npart) [60,61]: 



with Upp being the charged particle multiplicity in non-single diffractive pp interactions at 
the respective ^/s. Using parameters extracted from a Glauber calculation and fitting the 
parameter x to the PHOBOS data this analysis yields a fraction of 37% of the produced 
particles stemming from hard processes at =130 GeV. The centrality dependence of 
the charged particle multiplicity is nicely reproduced as well - see the left frame figure 5 
(taken from [59]). However, an analogous analysis utiUzing the saturation approach of 
the Color Glass Condensate yields similar agreement with the data, when employing a 
centrahty dependent saturation scale Qs{b) [60]: 



Since both, the conventional QCD and the high-density QCD (i.e. saturation) approach 
are able the describe the data, it is currently difficult to distinguish between these two 
approaches. Comparisons to other data, e.g. transverse momentum spectra, will be able to 
resolve this ambiguity. 




impart) 
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Figure 5. Left: Charged particle multiplicity per participant as a function of centrality 
measured by the PHOBOS experiment compared to various theoretical calculations. 
Right: tt" transverse momentum spectrum measured by PHENIX compared to pQCD 
based calculations employing different scenarios of partonic energy loss. 

Partons created in a hard scattering can act as very efficient probes of a QGP: a parton 
traversing a color confined medium of hadrons sees a relatively transparent system. How- 
ever, a parton passing through a hot medium of deconfined degrees of freedom is expected 
to loose a large amount of energy via gluon radiation. Since the radiated gluons may cou- 
ple to each other, the energy loss is predicted to be proportional to the square of the length 
traversed in the deconfined medium [63] - this phenomenon is commonly being refeiTed 
to as jet-quenching. RHIC experiments have established a significant suppression of high- 
Pt hadrons produced in central Ah-A collisions compared to those produced in peripheral 
Ah-A or binary scaled pH-p reactions, indicating a strong nuclear medium effect, compati- 
ble with the jet-quenching predictions [64,65]. Early calculations predicting a fairly large 
rate of energy loss around 2 GeV/fm have been shown to be incompatible with the mea- 
sured spectra at RHIC [64,65], whereas newer calculations based on the GLV formalism 
[66] for "thin" plasmas and including finite kinematic effects provide a far better agree- 
ment with the data [67] - similar to the assumption of an rather small energy loss around 
0.25 GeV/fm. 

However, alternative scenarios for partonic energy loss have to be considered: the parton 
may change its momentum through hard scattering in the deconfined medium [68] or after 
hadronizing through soft hadronic rescattering [69]. Figure 5 (taken from [68]) shows a 
parton cascade calculation which compares standard parton energy loss calculations with 
a collisional energy loss scenario - both approaches are compatible with the PHENIX data 
[64]. 



5.3 Balance Functions and Fluctuations 

While there is no doubt about the importance of partonic degrees of freedom for the ini- 
tial, early, reaction stages, one of the crucial questions is how long the deconfined state 
may have actually existed and whether this time-span is long enough for thermalization 
and collective effects to occur Balance functions offer a unique model-independent for- 
malism to probe the time-scales of a deconfined phase and subsequent hadronization [70]. 
If a long-lived QGP has been formed, a large number of quark-antiquark pairs need to be 
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Figure 6. Left: jK balance functions calculated for different hadronization 
times and scenarios. Right: Transport model calculation of charged particle fluctua- 
tions compared to analytical predictions for a resonance gas and a QGP. 



created close to the hadronization time, mainly due to entropy conservation constraints. 
This late-stage production of quarks could be attributed to three mechanisms: formation of 
hadrons from gluons, conversion of the non-perturbative vacuum energy into particles, or 
hadronization of a quark gas at constant temperature. Hadronization of a quark gas should 
approximately conserve the net number of particles due to the constraint of entropy conser- 
vation. Since hadrons are formed of two or more quarks, creation of quark-antiquark pairs 
should accompany hadronization. All three mechanisms for late-stage quark production 
involve a change in the degrees of freedom. Therefore, any signal that pinpoints the time 
where quarks first appear in a collision would provide valuable insight into understanding 
whether a novel state of matter has been formed and persisted for a substantial time. 

The link between balance functions and the time at which quarks are created has a simple 
physical explanation. Charge-anticharge pairs are created at the same location in space- 
time, and are correlated in rapidity due to the strong collective expansion inherent to a 
relativistic heavy ion collision. Pairs created earlier can separate further in rapidity due 
to the higher initial temperature and due to the diffusive interactions with other particles. 
The balance function, which describes the momentum of the accompanying antiparticle, 
quantifies this correlation. The balance function describes the conditional probability that 
a particle in the bin pi will be accompanied by a particle of opposite charge in the bin p^- 

B{p2\pi) = ^ {p{b,p2\a,pi) - p{b,p2\b,pi) + p{a,p2\b,pi) - p{a,p2\a,pi)} , 

(3) 

where p{b,p2\a,pi) is the conditional probability of observing a particle of type b in binp2 
given the existence of a particle of type a in bin pi . The label a might refer to all negative 
kaons with b referring to all positive kaons, or a might refer to all hadrons with a strange 
quark while b refers to all hadrons with an antistrange quark. 

The left frame of figure 6 shows /K^ balance functions as predicted in a simple 
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Bjorken thermal model for two hadronization temperatures, 225 MeV and 165 MeV as well 
as for fragmenting strings (utilizing PYTHIA), which would represent the hadronization 
scenario of a hadronic/string picture similar to that of RQMD or UrQMD. Since particles 
from cooler systems have smaller thermal velocities, they are more strongly correlated in 
rapidity and result in narrower balance functions. A strong sensitivity to the hadronization 
temperature and time can be clearly observed. The STAR collaboration has recently pub- 
lished first results on a charged particle balance function analysis at RHIC [71]: a striking 
centrality dependence is observed, with the width of the balance function decreasing as 
a function of increasing centrality, in line with what to expect in the case of a long-lived 
deconfined phase. 

The size of the average fluctuations of net baryon number and electric charge in a finite 
volume of hadronic matter differs widely between the confined and deconfined phases, 
due to the different value of the elementary charge carried by quarks vs. hadrons. These 
differences may be exploited as indicators of the formation of a quark-gluon plasma in 
relativistic heavy-ion colUsions, if the fluctuations created in the initial state survive un- 
til freeze-out [73,74]. The right frame of figure 6 (taken from [75]) compares transport 
calculations in the framework of a string/hadron model for SPS and RHIC with analyti- 
cal predictions for charge fluctuations in a hadronic resonance gas and a parton gas [75]. 
Since the transport model does not contain any deconfined degrees of freedom, the numer- 
ical calculation agrees well with the resonance gas prediction. 

Interestingly, sofar all experimental analysis for SPS and RHIC agree with the hadron 
gas prediction, giving no indication at all about a possible deconfined phase. Since many 
other observables are compatible with the assumption of deconfinement, it may very well 
be that the dynamics of hadronization strongly affects the charge fluctuation observables 
and masks the deconfined phase. 

Under certain assumptions, charge fluctuations can be directly expressed in terms of 
balance functions [72]: 



The size of the correction term is usually less than 5 % for electrical charge fluctuations, 
since the number of produced charges is much larger than the net charge. However, for 
baryon number fluctuations the additional term is not negligible, even at RHIC. 

Comparing the apparent lack of QGP indications of the charge fluctuation observables 
with the promising results of the balance function analysis one is lead to speculate whether 
relevant dynamical information of the two particle correlation is averaged out by integrat- 
ing over the correlation function in order to obtain the fluctuation observable. 

5.4 Collective Flow 

The study of deconfinement and a subsequent phase-transition to deconfined hadronic mat- 
ter poses a great challenge to microscopic transport models. In most approaches hadroniza- 
tion is a uni-directional process and the equation of state of the system ill-defined. One 
possible remedy is to use a hydrodynamical approach for the early deconfined phase of the 
reaction and subsequent phase-transition, coupled with a microscopic calculation for the 
later, hadronic, reaction phase in which the hydrodynamical assumptions are not vaUd any 
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Figure 7. Left: Inverse hadron rriT slopes at yc.m. = 0. The lines depict pure hy- 
drodynamics whereas the symbols refer to data and hydro+micro calculations. Right: 
collision number distributions for different baryon species in the hadronic phase of a 
hydro+micro calculation. 



longer [21]. In the following, such a combined macro+micro model will be used to study 
elliptic flow as well as the flavor- and mass-dependence of hadronic slope parameters (ie. 
radial flow): 

Radial and elliptic flow in non-central heavy ion collisions can constrain the effective 
Equation of State(EoS) of the excited nuclear matter. It has been shown that for an EoS 
with a first order phase transition, the above mentioned hybrid macron- micro models repro- 
duce both the radial and elliptic flow data at the SPS [21,76]. The centrality dependence 
of the elliptic flow coefficient V2 exhibits a strong sensitivity to the underlying EoS, which 
may help to constrain the EoS [77]. In addition, a number of features of the RHIC data can 
also be explained within the hybrid approach [77]: the observed elliptic flow and its de- 
pendence on pt and mass, the anomalous p/tt" ratio for pt ~ 2.0 GeV, and the difference 
in the impact parameter dependence of the p and (6 slope parameters. For an EoS without 
the hard and soft features of the QCD phase transition, the broad consistency with the data 
is lost. 

The left frame of figure 7 displays the inverse slope parameters T obtained by an ex- 
ponential fit to dNi/ cPrriTdy in the range rriT — rrii < 1 GeV for SPS and RHIC in a 
hybrid hydroH- micro calculation and compares them to SPS data [31]. The trend of the 
data, namely the "softer" spectra of E!'s and Q's as compared to a linear T(m) relation 
is reproduced reasonably well. This is in contrast to "pure" hydrodynamics with kinetic 
freeze-out on a common hypersurface (e.g. the T = 130 MeV isotherm), where the stiff- 
ness of the spectra increases linearly with mass as denoted by the lines in fig. 7. When 
going from SPS to RHIC energy, the model discussed here generally yields only a slight 
increase of the inverse slopes, although the specific entropy is larger by a factor of 4-5 ! 
The reason for this behavior is the first-order phase transition that softens the transverse 
expansion considerably. 

The softening of the spectra is caused by the hadron gas emerging from the hadroniza- 
tion of the QGP being almost "transparent" for the multiple strange baryons. This can 
be seen by calculating the average number of coUisions different hadron species suffer in 
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Figure 8. Left: Freeze-out transverse radius and time distribution dP N / {rTdrrdt fr) 
for pions at RHIC. Right: 7?out /Rside for RHIC initial conditions, as a function of Kt 
at freeze-out (symbols) and at hadronization (Unes). 



the hadronic phase (see right frame of fig. 7): whereas fi's suffer on average only one 
hadronic interaction, A'''s and A's suffer approx. 5-6 collisions with other hadrons before 
they free/e-out. 

Thus, one may conclude that the spectra of S's and especially fi's are practically unaf- 
fected by the hadronic reaction stage and closely resemble those on the phase boundary. 
They therefore act as probes of the QGP expansion prior to hadronization and can be used 
to measure the expansion rate of the deconfined phase. 



5.5 Two Particle Interferometry 

Kinetic freeze-out, at which all (elastic) interactions cease, terminates the evolution of the 
reaction discussed in figure 2. However, the freeze-out time of the system is ill-defined, 
since interactions cease locally on a particle per particle basis. The left frame of fig. 8 
shows the distribution of freeze-out points of pions in the forward light-cone, as well as the 
hadronization hypersurface from where all hadrons emerge in a hydroH-micro approach. 
Freeze-out is seen to occur in a four-dimensional region within the forward light-cone 
rather than on a three-dimensional "hypersurface" in space-time. Similar results have also 
been obtained within other microscopic transport models [78] when the initial state was 
not a QGP. It is clear that the hadronic system disintegrates slowly (as compared to e.g. the 
hadronization time), rather than emitting a "flash" of hadrons (predominantly pions) in an 
instantaneous decay. 

A first order phase transition leads to a prolonged hadronization time as compared to a 
cross-over or ideal hadron gas with no phase transition, and has been related to unusually 
large Hanbury-Brown-Twiss (HBT) radii [79-81]. The phase of coexisting hadrons and 
QGP reduces the "explosivity" of the high-density matter before hadronization, extending 
the emission duration of pions [79-81]. This phenomenon should then depend on the 
hadronization (critical) temperature Tc and the latent heat of the transition. For recent 
reviews on this topic we refer to [82,83]. 
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Figure 9. Left: Invariant mass spectrum of low mass dileptons calculated with and 
without hadronic medium modifications, including possible contributions from thermal 
dileptons stemming from a QGP. Right: Mass spectrum of intermediate mass dileptons 
- around M^e — 2.5 GeV the contribution of QGP radiation dominates above the 
hadronic background and Drell-Yan processes. 

It has been suggested that the ratio i?out / ^sido should increase strongly once the initial 
entropy density Si becomes substantially larger than that of the hadronic gas at [81]. 
The strong Tc dependence of i?out / ^sidc in such a purely hydrodynamical scenario can 
be seen in the right frame of fig. 8 (solid and dotted lines). Here, the i?out/J^sido ratio at 
hadronization (calculated in a hydrodynamical scenario) is compared to the ratio after sub- 
sequent hadronic rescattering and freeze-out in a combined hydroH-micro approach [84]. 
Clearly, up to Kt ^ 200 MeV (with Kt ~ {pi,t + P2,t)/2) i?out/^sidc is independent 
of Tc, if hadronic rescatterings are taken into account. Moreover, at higher Kt the de- 
pendence on Tc is even reversed: for high Tc the i?out / ^sidc ratio even exceeds that for 
low Tc- A higher Tc speeds up hadronization but on the other hand prolongs the dissipa- 
tive hadronic phase that dominates the HBT radii. This is because during the non-ideal 
hadronic expansion the scale of spatial homogeneity of the pion density distribution in- 
creases, as the pions fly away from the center, but the transverse flow can hardly increase 
to counteract. Therefore, after hadronic rescattering i?out / i?sidc does not drop towards 
higher Kt (in the range Kt < Sm^r). 

For central collisions of Au nuclei at ^/s = 130A GeV, data from STAR gives 
-Rout/^sidc — 1.1 at small Kt [85]. With increasing Kt the data remain flat or even 
decrease slightly - a trend which is currently not understood theoretically and cannot be 
reproduced by any of the dynamical models on the market. Kaon interferometry may pro- 
vide crucial insight into this so-called HBT-puzzle at RHIC [86]: the kaon phase-space 
density is much lower than that of the pions, thus dramatically reducing higher order cor- 
relation effects. In addition, high-pf kaons are emitted to a larger fraction directly from the 
phase-boundary, increasing the sensitivity to the QGP equation of state. 



5.6 Dileptons and Charm 

First dilepton and charm measurements will soon be forthcoming at RHIC and are highly 
anticipated: in the dilepton sector, current theoretical expectations are that in the low mass 
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region thermal dileptons from a QGP will be strongly suppressed compared to hadronic 
contributions - mainly from in-medium pion-pion annihilation. Current state of the art 
calculations point towards a strong broadening of vector mesons in medium [87] - the p and 
a; peaks should be almost completely dissolved and the significantly broadened (see left 
frame of figure 9, taken from [87]). At intermediate invariant masses (around 2 - 2.5 GeV) 
thermal dileptons radiated from a QGP start to outshine hadronic sources [88,89,87], while 
for higher masses the yield of Drell-Yan processes exceeds that of thermal dileptons from 
a QGP (see right frame of figure 9). Open charm may give a large contribution to the 
spectrum in the intermediate mass range as well - however independent measurements of 
open charm will allow to subtract that contribution from the invariant mass spectrum and 
thus will allow a measurement of the thermal contribution [90,87]. 

In the charmonium sector, predictions on what to expect at RHIC vary widely - from 
total suppression to a strong enhancement of the J/i/; and yield! Originally it was pre- 
dicted that the suppression of heavy quarkonia-mesons could provide one of the signatures 
for deconfinement in QCD at high temperatures [91]. The idea was based on an analogy 
with the well known Mott transition in condensed matter systems. At high densities, Debye 
screening in a quark-gluon plasma reduces the range of the attractive force between heavy 
quarks and antiquarks, and above some critical density screening prevents the formation of 
bound states. The larger bound states are expected to dissolve before the smaller ones as 
the temperature of the system increases. The il/ and Xc states are thus expected to become 
unbound just above T^, while the smaller ^ state may only dissolve above « 1.2Tc. Heav- 
ier bb states offer the same features as cc states, but require much shorter screening lengths 
to dissolve [92]. The T{bh) state may dissolve only around 2.5 T^, while the larger excited 
T' could also dissolve near Tc. 

At SPS energies, less than one charmonium is on average produced per heavy-ion reac- 
tion [93,94] - here the screening argument is most compelhng and the debate has shifted 
towards competing hadronic processes [95]: the size of the hadronic absorption cross sec- 
tions for charmonia remains a contentious issue. However, these cross sections are nec- 
essary to estimate the contribution of the hadronic phase to charmonium suppression and 
unless a more solid understanding of these cross sections is obtained, the interpretation of 
the SPS charmonium data will remain under debate. At RHIC, several cc pairs may be 
created in one event - thus allowing for a coalescence of a cc pair which was not produced 
in the same hard process. Estimates taking this effect into account actually predict an en- 
hancement of charmonium production in a QGP, compared to a purely hadronic scenario 
[96]. 

6. Outlook 

Thanks to the vast amount of exciting new RHIC data, QGP theory is currently the most 
active, dynamic and exciting field in theoretical nuclear physics. Already the current data 
are of such high quality that they might allow not only the observation of a deconfined 
phase of matter, but also the characterization of the QGP equation of state! In order to 
accomplish this feat, a consistent picture of all phases of a heavy-ion reaction needs to be 
developed, not limited to a few chosen observables. Proper baselines for comparison need 
to be estabhshed as well - the importance of the RHIC p-np and p-nA program cannot be 
overstated: these experiments will serve to constrain theories and models as well as help 
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to generate a better understanding of initial state effects. 

There are a number of milestones which QGP theory needs to address in the forthcoming 
years to fulfill the promise the data is currently showing us: 

• the development of reliable lattice gauge calculations at finite temparature and chem- 
ical potential 

• the determination of hadron properties at high densities and temperatures 

• generating a better understanding on the mechanisms of hadronization, both, for 
elementary hadron-hadron reactions as well as for bulk QCD matter 

• the Unk-up of initial state saturation models to the later-stage dynamics of a heavy- 
ion colUsion, including pQCD processes 

• the investigation of mechanisms for chemical equilibration, both in the confined as 
well as in the deconfined phase 

• the development of a transport theory based on QCD which treats both, hard as well 
as soft processes consistently. 
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